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Abstract 



The polarizations of the Y(1S), Y(2S), and Y(3S) mesons are measured in proton- 
proton collisions at \/s = 7 TeV, using a data sample of Y(nS) — > u + u ~ decays col- 
lected by the CMS experiment, corresponding to an integrated luminosity of 4.9 fb _1 . 
The dimuon decay angular distributions are analyzed in three different polarization 
frames. The polarization parameters A#, A^, and A^g,, as well as the frame-invariant 
quantity A, are presented as a function of the Y(nS) transverse momentum between 
10 and 50GeV, in the rapidity ranges \y\ < 0.6 and 0.6 < \y\ < 1.2. No evidence of 
large transverse or longitudinal polarizations has been seen in the explored kinematic 
region. 
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Studies of heavy-quarkonium production play an important role in the detailed investigation 
of quantum chromodynamics (QCD). Given their high mass, heavy-quarkonium states are ap- 
proximately nonrelativistic systems, allowing the application of theoretical tools that simplify 
and constrain the analyses of nonperturbative effects [1]. The differential cross sections of ]/ip 
and Y mesons produced at Tevatron [2-4] and LHC [5-7] energies can be reproduced by cal- 
culations based on nonrelativistic QCD (NRQCD) [8|, dominated by "color octet" production. 
However, the corresponding predictions [9| of strong transverse polarizations (dominant an- 
gular momentum component J z = ±1 with respect to the quarkonium momentum direction) 
are in stark disagreement with the negligible polarizations measured for the J/t/> |flO|. The Y 
satisfies the nonrelativistic approximation much better than the ]/xp, making the Y polarization 
a more decisive test of NRQCD, especially at asymptotically large transverse momentum, pj. 
The existing measurements, however, are inconclusive, with the CDF [11] and DO lTl2~| results 
in mutual contradiction. 

The polarization of the (J PC = 1 ) Y states can be measured through the study of the angular 
distribution of the leptons produced in the Y — > decay [13 1, 

W(cos #, <»| A) oc — — (1 + A#cos 2 # + A™ sin 2 #cos2<» + sin2#cos a) , (1) 

(3 + A#J 

where i? and cp are the polar and azimuthal angles, respectively, of the ]d + with respect to the 
z axis of the chosen polarization frame. As pointed out in Refs. ||T3T[T7| , improved experimen- 
tal measurements of quarkonium polarization require measuring all the angular distribution 
parameters, A = (A#, A^, A#«,), in different polarization frames, as well as a frame-invariant 
polarization parameter, A = (A# + 3 (1 — A^). This approach has already been followed in 
the Y(nS) polarization analysis of CDF [11], and in some recent theory calculations [fl8| . 

This Letter presents the first measurement of the polarizations of the Y(nS) mesons produced 
in pp collisions at a center-of-mass energy of 7 TeV. The analysis is based on a dimuon sample 
collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 
4.9 fb" 1 and containing 252 000 Y(1S), 94 000 Y(2S), and 58 000 Y(3S) mesons (after all selection 
criteria). 

The analysis uses an unbinned likelihood approach, completely independent of theoretical hy- 
potheses on the production kinematics, and with the lepton efficiencies as the only external 
input. The results are obtained in three frames, with different directions of the quantization 
axis: the center-of-mass helicity (HX) frame, where the polar axis coincides with the direction 
of the Y momentum; the Collins-Soper (CS) frame |Tl9| , whose axis is the average of the two 
beam directions in the Y rest frame; and the perpendicular helicity (PX) frame [20 1, orthogonal 
to the CS frame. The y axis of the polarization frame is taken, in all cases, to be in the direction 
of the vector product of the two beam directions, Pi x P? and P2 x Pi for positive and negative 
rapidity, respectively. 

The central feature of the CMS apparatus |2l"1[ is a superconducting solenoid of 6 m internal 
diameter, providing a 3.8 T field. The main subdetectors used in this analysis are the silicon 
tracker and the muon system. The silicon tracker, composed of pixel and strip detector mod- 
ules, is immersed in the magnetic field and enables the measurement of charged-particle mo- 
menta over the pseudorapidity range \tj\ < 2.5. Muons are measured in the range \r/\ < 2.4 
using gas-ionization detectors embedded in the steel return yoke of the magnet and made 
using three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. 
The events were collected using a two-level trigger system. The first level consists of custom 
hardware processors and uses information from the muon system to select events with two 
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muons. The "high-level trigger" requires an opposite-sign muon pair with invariant mass 
8.5 < M < 11.5 GeV, |y| < 1.25, pj > 5 or 7GeV (depending on the instantaneous luminosity), 
and vertex fit % 2 probability greater than 0.5%. 

In the offline analysis, dimuons are formed by combining pairs of opposite-sign muons (tracks 
in the silicon tracker matched to tracks in the muon detectors) that satisfy several quality cri- 
teria, including the number of tracker hits, the muon-track fit quality, and the vicinity of the 
track to the closest primary vertex along the beam line. The selected muons are required to 
satisfy \t]\ < 1.6 and to have pj above 4.5, 3.5, and 3.0GeV for \rj\ < 1.2, 1.2 < \t]\ < 1.4, and 
1.4 < \t]\ < 1.6, respectively, to ensure accurately measured muon detection efficiencies. Sub- 
sequent to the offline trigger confirmation, the combinatorial background from uncorrelated 
muons is reduced by requiring a dimuon vertex fit x 2 probability larger than 1.0% and a dis- 
tance between the dimuon vertex and the closest primary vertex less than twice its uncertainty. 
The analysis is performed in five dimuon pj bins, of edges 10, 12, 16, 20, 30, and 50 GeV, and 
two |y| ranges, 0.0-0.6 and 0.6-1.2. 
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Figure 1: Dimuon mass distributions in the Y region for \y\ < 0.6 (open squares) and 
0.6 < |y| < 1.2 (closed circles). 

The dimuon mass distribution, shown in Fig. [lj is well described by three Crystal-Ball func- 
tions |22] representing the Y peaks, and by a second-order polynomial function determined 
from the low- and high-mass sidebands, located below the Y(1S) and above the Y(3S), respec- 
tively. The dimuon mass resolution is better than 70MeV for \y\ < 0.6, increasing to 95MeV 
in the 0.6 < \y\ < 1.2 range, where the Y(2S) and Y(3S) peaks partially overlap. Within a 
±1 standard deviation (o~) window around the Y(nS) masses, the cross-feed between the Y(2S) 
and Y(3S) is below 4%, and the background fractions are 4-8%, 9-18%, and 12-28% (increasing 
with decreasing pj), for the Y(1S), Y(2S), and Y(3S), respectively. 

The single-muon detection efficiencies are measured with a "tag-and-probe" technique |23"| 
using event samples collected with dedicated triggers enriched in dimuons from ]/xp decays. 
The trigger and reconstruction efficiencies must be accurately determined to avoid biases on 
the angular distributions, which could mimic polarization effects. The technique has been 
validated in the fiducial region of the analysis with detailed Monte Carlo (MC) simulation 
studies. The single-muon efficiencies are measured and parametrized as a function of pj in 
eight tj bins. Their uncertainties, reflecting the statistical precision of the calibration samples 
and possible imperfections of the parametrization, contribute to the systematic uncertainty on 
the final results. The dimuon trigger and the selection criteria applied at the dimuon level 
could potentially introduce muon-pair correlations, making the dimuon detection efficiencies 
different from the product of the efficiencies of the two single muons. Detailed MC simulations 
show that such correlations are essentially independent of cos # and <p, in the phase space 
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selected for the measurement. Residual effects are incorporated into the systematic uncertainty. 

A fit to the dimuon mass distribution provides the fraction of background events, /b, under 
each of the three Y mass peaks, for a given definition of the signal region. The angular distribu- 
tions of these background events are modeled as weighted sums of the distributions measured 
in the mass sidebands (defined with negligible signal contamination), with weights derived 
under the assumption that they change monotonically with dimuon mass. The background 
dimuons are subtracted on an event-by-event basis using the likelihood ratio Lb/ L(s+bV where 
Lb is the likelihood of an event to be background, reflecting the background model, and L( S+B ) 
is its likelihood to be either signal or background, reflecting the data distribution. Both likeli- 
hoods are functions of the variables pj, |i/|,M,cos#, and (p. They have relative normalization 
equal to /b, so that a fraction /b of events distributed according to the (pj, \y\,M, cos d, q>) dis- 
tribution of the background model is removed from the data sample. 

The posterior probability distribution (PPD) for the average values of the Y polarization pa- 
rameters (A) inside a particular kinematic cell is then defined as a product over the remaining 
("signal-like") events (z), 

H*) = Yl £ (P? >Pi ] ) > ( 2 ) 

i 

where £ represents the event probability distribution as a function of the muon momenta pi 2 
in event i. Unlike most polarization analyses, we do not use simulated (cos q>) acceptance 
and efficiency maps, averaged over all events in the considered kinematic cell. Instead, the 
procedure exploits the efficiency measurement as a function of muon momenta, attributing to 
each event a probability dependent on the full event kinematics (not only cos 1? and cp) and on 
the values of the polarization parameters. The event probability is defined as 

£ (Pi, pi) = W(cos #, <p|A) e(pi, p 2 ) , (3) 

W {A) 

where e(p\, pi) is the detection efficiency. The normalization factor A/"(A) is calculated by inte- 
grating W • e over cos 1? and q> uniformly, using (pj, \y\, M) distributions determined from the 
background-subtracted data. 

The background subtraction procedure is repeated 50 times to evaluate the statistical fluctua- 
tions associated with its random nature, and the PPD is obtained as the average of the PPDs 
calculated over the 50 ensembles of signal-like events. 

The analysis framework, including the effects of the detection efficiencies, has been tested with 
pseudo-experiments based on simulated samples. Each test involves 50 pseudo-experiments 
and evaluates a specific systematic uncertainty. The pseudo-samples are individually gener- 
ated and reconstructed, leading to statistically independent determinations of the polarization 
parameters. The difference between the median of the 50 results and the injected polariza- 
tion parameters provides the systematic uncertainty corresponding to the effect under study. 
The precision of the method to extract the signal polarization is evaluated for several signal 
and background polarization scenarios. The influence of a possible residual bias from muon 
or dimuon efficiencies, stemming from the tag-and-probe measurement precision or from the 
efficiency parametrization, is evaluated by applying uncertainty-based changes to the efficien- 
cies used in the extraction of the polarization parameters. The monotonicity hypothesis in the 
data-driven modeling of the background angular distribution under the Y(nS) peaks has been 
tested by varying the signal region from ±3c to ±1 a around the Y(nS) masses (with corre- 
sponding kinematic corrections determined from simulation). Despite significant changes in 
/b (from 40% to 28% for the Y(3S) at low pj and |y| < 0.6, for instance), the results remain 
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Figure 2: Correlations between A<p and A# (a), and A#,p and A<p (b) for Y(1S) with |y| < 0.6 and 
30 < p T < 50GeV. The 68.3% and 99.7% CL contours are shown for the CS and HX frames. 
The shaded areas represent physically forbidden phase space regions |17|. 

essentially identical for all three states. Larger variations are observed by modifying the rel- 
ative weights of the low- and high-mass sidebands in the background model composition. A 
conservative range of hypotheses is considered, such as assuming that the background under 
the Y(1S) (Y(3S)) peak resembles exclusively the low-mass (high-mass) sideband, or assuming 
that it is reproduced by an equal mixture of the two sideband distributions. While there is no 
dominant source of systematic uncertainty in the Y(1S) case, the total systematic uncertainty 
of the Y(2S) and Y(3S) states is dominated by the background model uncertainty, especially at 
low pj. At high pj, the statistical uncertainties dominate. For example, the statistical uncer- 
tainties in Atf(PX) at \y\ < 0.6 for the Y(1S) (Y(3S)) are of order 0.1 (0.2) at both low and high 
pi; the corresponding systematic uncertainties have a similar magnitude at low pj and are a 
factor of two (three) smaller at high pj. 

Each final PPD is an envelope of all PPDs corresponding to the ranges of hypotheses covered 
by the systematic uncertainties. One- and two-dimensional projections of each final PPD are 
calculated by numerical integration. The highest posterior probability in each one-dimensional 
projection is used to estimate the best value of the associated polarization parameter. Inter- 
vals [Ai, A2] corresponding to a given confidence level (CL), are calculated by identifying two 
regions of the parameter space, [— oo,Ai] and [A2,oo], each containing 0.5 • (1 — CL)% of the 
one-dimensional projection of the PPD. Figure [2] shows two projections of the final PPD for the 
Y(1S) at |y| < 0.6 and 30 < p T < 50GeV, displaying the 68.3% and 99.7% CL contours for the 
CS and HX frames. 

Figure [3] shows, for the rapidity range 0.0-0.6, one-dimensional profiles (68.3%, 95.5%, and 
99.7% CL intervals) of the PPDs of the parameters A#, A<p, and A^, for the Y(1S), Y(2S), and 
Y(3S) states, in the HX frame. Similar values are obtained in the 0.6-1.2 rapidity range [24J. 
Figure [i] displays the corresponding results for the frame-invariant parameter A, including also 
the CS and PX values. The results obtained in the three frames are in good agreement, as 
required in the absence of unaccounted for systematic effects. Complete tables of results for 
A#, A,p, Afig,, and A, for the three Y states and in the three frames considered in this analysis, are 
available in Ref. ||24| . 

All the polarization parameters are compatible with zero or small values in the three polar- 
ization frames, excluding that a significant polarization could remain undetected because of 
smearing effects induced by unfortunate frame choices. The indication that the Y(nS) reso- 
nances are produced as an unpolarized mixture might be related to the fact that the measure- 



5 



1 




Y(1S) ; 


Y(2S) i 




Y(3S) 


0.5 

K ° 


MM-M--- m ■ _ 




■ ■ ■ ■ ■ 

■ ■ ■ 


■ III 1 













-0.5 
-1 


- CMS pp Vs = 7TeV L 


= 49 f b" 1 


uy frame li/l fl fi - 
riA name, \y\ <. u.o 






0.4 
0.2 

\ 


4 Stat, uncert., 68.3 % CL 
Tot. uncert., 68.3 % CL 
Tot. uncert., 95.5 % CL 
Tot. uncert., 99.7 % CL 


Y(1S) : 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 
Y(2S) : 


1 M II | 1 ! II | II II | II II | 1 II ! | II 


1 1 | 1 1 1 1 | 1 1 1 1 

Y(3S) 






MM - m --- jg 


■ 




-0.2 












-0.4 


I " " I " " I " ' ' I " " I ' i " I " " 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 


0.4 




Y(1S) : 


Y(2S) : 




Y(3S) 


0.2 














■)-■'■■»-, 1-- 




1 


H-I---V 1 






t | I 




-0.2 










-0.4 













10 15 20 25 30 35 40 45 10 15 20 25 30 35 40 45 10 15 20 25 30 35 40 45 



P T [GeV] P T [GeV] P T [GeV] 

Figure 3: Values of the A# (top), Aq, (middle), and A#«, (bottom) parameters for the Y(1S) (left), 
Y(2S) (middle), and Y(3S) (right), in the HX frame, as a function of the Y p T for \y\ < 0.6. The 
error bars indicate the 68.3% CL interval when neglecting the systematic uncertainties. The 
three bands represent the 68.3%, 95.5%, and 99.7% CL intervals of the total uncertainties. 
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Figure 4: Values of A for the Y(1S), Y(2S), and Y(3S) states (left to right), in the HX, CS, and 
PX frames, for the \y\ < 0.6 (left) and 0.6 < \y\ < 1.2 (right) ranges. The bands and error bars 
have the same meaning as in the previous figure. 

ments do not distinguish directly produced Y mesons from those produced in the decays of 
heavier bottomonium states. 

In summary, the polarizations of the Y(nS) mesons produced in pp collisions at = 7TeV 
have been determined as a function of the Y pj in two rapidity ranges and in three different 
polarization frames, using both frame-dependent and frame-independent parameters. The re- 
sults exclude large transverse or longitudinal Y(nS) polarizations, beyond the pj and y ranges 
probed by previous experiments, especially for the Y(3S) state, and tend to be in disagreement 
with theoretical expectations for high-energy hadron collisions H9l l25| . 
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